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Abstract 

We properly define off-shell K — > tt transition amplitudes and use them 
to extract information for on-shell K — » tttt amplitudes within Chiral Per- 
turbation Theory. At order p 2 in the chiral expansion all three param- 
eters of weak interaction can be determined. At order p 4 we are able 
to fix eleven additional constants out of thirteen contributing to off-shell 
K —* tt transitions, which leaves four undetermined constants in the on- 
shell K — > tttt amplitudes. All 0(p ) contributions have been exactly 
derived with ^ 0. We finally discuss the weak mass term issue and 
find contributions to on-shell AS = ±1 Kaon decays, in particular to tran- 
sitions like Kl — ► 77, — > n + n~ and K$ — > 7r°77 at the lowest non-zero 
order. 
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1 Introduction 



The explanation of the A J = 1/2 rule in K — > tttt decays remains one of the chal- 
lenges in Kaon physics and in our understanding of strong interactions. Various 
non-leptonic Kaon decays are also used to put limits on CP-violation and several 
other quantities of the Standard Model and extensions of it ||. The short dis- 
tance part of the relevant operators can be treated using renormalization group 
within perturbative QCD, while the computation of matrix elements of the rele- 
vant operators between meson states is a pure non perturbative problem. 

In the long term lattice QCD should be able to perform a direct computation 
of weak matrix elements. It is however much easier on the lattice^], and often 
also in analytical attempts to reproduce the weak matrix elements, to calculate 
correlators involving fewer external legs. As a first step, current algebra can be 
used to relate K — > 2tt to K — > tt amplitudes where it involves off-shell K — > tt 
transitions. Chiral Perturbation Theory (CHPT) || is however the more modern 
tool to exploit the consequences of current algebra. At lowest order in the chiral 
expansion this problem was first worked out in [|J and subsequently the non- 
analytic parts of the loop contributions to K — > tttt and K — > tt were calculated 
in |J . It was also discussed in || and in the context of Wilson Fermions on the 
lattice in 0. In this paper we extend the previous study in two ways: 

1) We systematically go to next-to-leading order (i.e. order p 4 ) in CHPT and 

2) Instead of the vague notion of off-shell Kaon and pion fields, we use pseudo- 
scalar current correlators which are well defined quantities^. 

The use of this type of correlators to extract information on non-leptonic 
Kaon matrix elements is quite common in lattice studies (see for instance 0), 
though in those cases an on-shell extrapolation is usually performed. In addition, 
this extrapolation is done at lowest order p 2 in the chiral expansion, i.e. using 
pure current algebra relations. However, due to the large Kaon mass, one expects 
non-negligible higher order CHPT corrections to Kaon weak matrix elements and 
in general to the pseudo-scalar current correlators involved. The use of the off- 
shell behaviour of this type of correlators to obtain additional information on the 
relevant matrix element has been advocated in ||, where it was used to unravel 
the quark mass dependence of Bk, and in |§ to disentangle the structure of the 
electromagnetic mass differences. 

That chiral corrections are important in non-leptonic Kaon decays is already 
known since a long time W\ and has been fully worked out in CHPT by 



n], [12], [13| . Here we present results for the octet and the 27-plet contributions 
to K — > tttt transitions both at order p 2 and order p 4 and without neglecting 
m 2 n ./m 2 K suppressed contributions. (We have a small disagreement here with 

1 Computing the K — > tttt amplitudes directly is quite difficult because of the Maiani- Testa 
argument ||. 

2 We have verified that in the cases discussed here the use of other two-point correlators does 
not yield additional information. 
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respect to previous literature for the 27-plet [|I3"|.) While for the physical case 
neglecting m 2 /m 2 K terms is a reasonable approximation it will not be the case 
for foreseeable lattice calculations. We keep m = m u = rrid 7^ throughout the 
derivation. Another issue to be clarified is the weak mass term contribution. It is 
well known that the weak mass term, which gives rise to the tadpole contributions 
in the original formulation of the weak effective Lagrangian [|TTJ, does not enter 
the K — > mi on-shell matrix elements at order p 2 [14], 15| . It does however 



contribute in a well defined way to off-shell quantities at order p 2 and higher. For 
this reason we shall discuss the precise role of the weak mass term up to order p 4 
in Kaon transition amplitudes in Sect. |6|. 

It turns out that, while at order p 2 all the weak parameters can be determined 
from our two-point correlators^, this is no longer true at order p A . There are in 
total nineteen parameters entering the weak effective Lagrangian up to that order 
as discussed in Sect. pTT| . Of these, we can obtain fourteen using our procedure. 
Five of them can be obtained in more than one place, thus providing as many 
CHPT relations. 

Several relations are also implied between different two-point correlators so 
that the same quantities can also be used to check how well calculations within 
CHPT obey the chiral symmetry predictions to order p A . In particular, we can 
obtain several coefficients of terms which involve quark masses at higher order. In 
the purely strong sector these are the most difficult ones to predict from models 
and/or dispersive constraints. Determining some of them through our procedure 
will provide a good check on models used in this context (see e.g. 

We study in CHPT the pseudo-scalar current correlators 

H ij (q 2 ) = i J d*ze**<0|T(i* t (0)P J >)e irAS =«) |0) (1.1) 

in the presence of strong interactions. Above, a = ±1, ±2 stands for | AS\ = 1,2 
transitions and i,j are light quarks combinations corresponding to the octet of 
light pseudo-scalar mesons: 

t /^ _ 1 („■.„. J- J\ . r»r+/_\_j. r,K°, 



P w (x) = —i= [ui'YsU — di'jsd) ; P 77 (x) = dij 5 u; P (x) = sij^d; 
V2 v 



P K+ (x) = sij 5 u; P m (x) = -y= (uijsu + di^d — 2sij 5 s^ . (1.2) 

The effective action of weak interactions Tas^o describes strangeness changing 
processes in one and two units. Within the Standard Model it can be written as 
follows 

T AS =a = -C AS =a G F J d 4 y AS =a(y) , (1-3) 

where A s=a is a sum over the effective operators arising after integrating out 
the heavy bosons, i.e. W, Z, and the Higgs boson, and heavy fermions, i.e. top, 

3 These are called three-point correlators in lattice QCD because of the extra weak vertex. 
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bottom, and charm quark (see e.g. [M [T7|, for the actual derivation). The 
constant C^s=a collects Clebsch-Gordan factors and Gp is the Fermi constant. 



In (O), the first term in the expansion of exp[zFAs= a ] describes strangeness 
zero changing transitions, the second term describes strangeness one and two 
changing processes, while the third term includes (AS = ±1) 2 transitions. The 
| AS 1 1 = 2 case relevant to the factor which parameterizes the K° — K° mixing 
was already studied in ||; in Section 3, we will just repeat the relevant expressions 
for completeness. The plan of the paper is as follows. In Section 2 we construct 
the weak effective Lagrangian up to order p 4 and relevant to our analysis. The 
1/N C counting of the weak constants is also done in subsection 2.2. In Section 
3 the fully renormalized two-point current correlators up to order p 4 are derived 
for \AS\ = 0, 1, 2 cases. The non-analytic contributions to one loop are collected 
in Appendix [A[ In Section 4 the K — > ixix on-shell amplitudes are derived in 
CHPT up to order p 4 . In Appendix |B| are the non-analytic contributions to one 
loop. Section |5| is devoted to the connection between off-shell K — > 7r transition 
amplitudes and on-shell K —>■ 2tt amplitudes. Resonance saturation also for the 
27-plet sector is used here and derived in Appendix |C|. Finally, in Section |6| we 
clarify the role of the weak mass term in Kaon decays up to order p 4 and in 
Section [7| we state our conclusions. 



2 CHPT Lagrangian and l/A^- Discussion 

At lowest order in CHPT (i.e. 0(p 2 )) the strangeness changing interactions up 
to two units amongst the pseudo-Goldstone bosons and external scalar, pseudo- 
scalar, vector and axial-vector sources (neglecting virtual photon interactions) 
are described by the following effective Lagrangian: 



r {2) r <zi 
Z -A5=0 ' Z -A5=l t ^AS 



(2) 



,(2) 



The first term is the strong interaction Lagrangian 



,(2) 
'AS=0 



-f [tr(u%)+ tr( X+ )] 



(2.1) 



(2.2) 



where tr(A) is the flavour trace of A, 

w M = i«t(D M £7)«t =w t (2.3) 

and U = u 2 = exp(iv / 2$/-Fo) is the exponential representation incorporating the 
SU(3) matrix of the octet of light pseudo-scalar mesons 

72 Vq 

X " P - + K° . (2.4) 



$(rr) 



V2 



71 



V 



K~ 



V6 J 
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D^U denotes the covariant derivative on the U field 

Dpi! = d^U - + aJU + iU{v„ - a M ) , (2.5) 

where v^{x) and a M (x) are external SU(3) vector and axial-vector matrices (no 
singlet component will be included in the present analysis). The matrix x+ m 
Eq. ( p.2|) and its pseudoscalar counterpart x~ are defined as follows 

X± = u] XU ] ± u\ ] u , (2.6) 

where x = 2B (M. + s(x) +ip(x)), s(x) and p(x) are external scalar and pseudo- 
scalar SU(3) sources and M. = diag(m u , m d , m s ) is the light quarks mass matrix. 
The constant B is related to the vacuum expectation value of the scalar quark 
density 

(0\qq\0) \ q=uAs = -F 2 B (1 + 0{M)) . (2.7) 

With this normalization, F Q is the chiral limit value of the pion decay constant 
F-x ~ 92.4 MeV. In the absence of the U(1) J 4 anomaly (i.e. in the large N c 
limit) ||19|| , the U(3) singlet field rji becomes the ninth Goldstone boson which is 
incorporated in the $(x) field as 

In this work we limit ourselves to the octet symmetry case, meaning that we 
assume the singlet degree of freedom r/i as heavy and integrated out. This is 
enough for our purpose of showing how to relate off-shell K — > 7r and K — > rjs 
transitions to on-shell K — > tctt amplitudes. In octet symmetry, i.e. det(w) = 1 
and tr(w M ) = 0, the AS = ±1 contribution to the l.h.s. of ( |1.3| ) is given by 



/"( 2 ) (~l 774 

*-AS=l — U r 



G 8 tr(A 32 «X)+ G 8 tr (A 32 x+) 



+G 27 t^ kl tr (AjjM M ) tr(A w u' 



+ h.c, (2.9) 



where i, j, k, I — 1, 2, 3 correspond to the light flavour indices u, d, s and the tensor 
t ij ' kl has 

£21,13 t 13 ' 21 = - - t 22 ' 23 = t 23,22 = — - ■ 

3 ' 6 ' 

t 23,33 = t 33,23 = _I. t 23,ll =t ll,23 = I £ 1Q) 

6 3 
and zero otherwise. The matrix Ay is defined as 

Ay = uX^v) 

( X ij)ab = hadjb, (2.11) 
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with Si a the Kronecker delta acting on the SU(3) light flavour space. The constant 
C = Cas=i hi (|1.3| ) includes normalization factors and the Cabibbo-Kobayashi- 
Maskawa matrix elements 

C = -\^V ud V: s . (2.12) 

The couplings G% and G' 8 modulate octet operators under SU(3)l x SU(3)#, 
while G27 modulates a 27-plet operator. For on-shell K — ► tttt transitions at 
order p 2 one can set G' 8 = M, [T4|, [T5| . At order p 4 this question was studied in 
11 , ^0] ; the result is that one can always use a basis for the order p 4 counterterms 
where the effects of the weak mass term (G' s ) on the on-shell K —>■ ttit amplitudes 
are fully reabsorbed at this order. Clearly, the use of the shifted basis implies 
a redefinition of the order p 4 couplings in order to absorb the weak mass term 
contributions. This was not done in . Since G' 8 does always appear in off-shell 
K — > 7r transitions, we keep the unshifted basis in our analysis, where K — > txtx 
amplitudes explicitly contain order p 4 contributions proportional to G' 8 . 
The I AS 1 1 = 2 term in Eq. (|1.3|) can be written as follows 



4 2 5=2 = C AS=2 F 4 G 27 tr(A 32Ufl ) tr(A 32 <) + h.c. , (2.13) 



with 



C AS =2 = ~^(ml ml M 2 W , V C km) (2.14) 

and Tim 2 , m 2 , Myy, Vckm) being a known function of the heavy fermions and 
bosons masses, and the Cabibbo-Kobayashi-Maskawa matrix elements [ffj. 

The weak couplings G$, G' 8 and G27 in Eqs. (|2.9| ) and ( p.!3| ) are dimensionless 
and they are related to those used in [Tlj] as follows 

CF 4 G s = c 2 CF 4 G 27 = 3c 3] C F 4 G' 8 = c 5 . (2.15) 
2.1 The order p A 

At next-to-leading order in the chiral expansion (i.e. at order p 4 ) the complete 
list of counterterms in the strong interaction sector and in the octet symmetry 
case has been given by Gasser and Leutwyler in | 2T| . In the SU(3) flavour case 



and for on-shell Green's functions there appear ten counterterms denoted with 
Li, % — 1, . . . , 10, while in the off-shell case there are two extra contact terms Hi 
and H 2 involving external sources only. 

The complete basis of counterterms in the weak interaction sector at order p 4 
and describing transitions with strangeness changing in one and two units was 
first derived by Kambor, Missimer and Wyler in |TT| and in P2|. This basis 



was confirmed and reduced to various minimal sets by Esposito-Farese for the 
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31, and by Ecker, Kambor and Wyler for the octet 
±1 transitions we use a minimal set of operators 



octet and 27-plet operators 
operators |24[ . 

For the analysis of the AS 
which differs from the one in [23| in one octet operator, but has the advantage 
of producing shorter expressions. Our octet subset coincides with that of Ecker 
et al. in [24] up to two operators. We also give below the translation from the 



counterterms we are using to those in |24]. 

In the octet symmetry minimal set of counterterms contributing to 

the K — > 7r and K — > tttt transitions at order p 4 in CHPT is given by: 



(4) 

AS=1 



E x O\ + E 2 G 8 + E 3 0\ + E 4 8 + E 5 8 



+ E 11 8 u + E 12 8 2 + E 13 O s 13 + E 15 



= C Eq g 8 

+ E 10 O 8 

+ C Fq G 2 7 

+ D 6 Of + D 7 Of +h.c. 



15 



D x Of + D 2 Of + D A Of + D 5 Of 



The octet operators above are 

Of = 



o\ 
o\ 
o\ 

Oto 

o 8 

^12 

o 8 

^13 

o 8 



The 27-plet operators are 

Of = 

Of = 

of = 

of = 

of = 

of = 



t ij,k 
t ij,k 
t ij,k 
t ij,k 
t ij,k 
t ij,k 



tr 
tr 
tr 
tr 
tr 
tr 
tr 
tr 
tr 
tr 

tr 
tr 
tr 
tr 
tr 
tr 



A32X+) tr(x+) ; 

A 32 x-x-) ; 

A32X-) tr(x-) ; 

a 32 bc+,x-]) ; 

A 32 { X +,u^ Ufl }) ; 

A 32 u^x+ u ^) ; 

A 32UfM ) tr(K,x+}) 
A32X+) tr ; 
A 32 • 



Kl 



Kl 



AijU^) 
AijU^) 



A 



Kl 



x+. 



A^) 



tr (A k ix+) ; 
tr (A klX -) ; 
tr {A kl K,x+}) 
tr(A H K,x-]) ; 

tr (Amu^) ; 
tr(A w O tr(x+ 



where the tensor t 1 -*' is defined in ( |2.10|) . 



(2.16) 



(2.17) 



(2.18) 



The translation from the octet counterterms in Q2.17 ) to the ones used by 



Ecker et al. \*M\ is as follows 

N 5 = E 10 



En 



N e = E n + 2E 



12 



6 



iV 7 


= 


; N 8 


= E n ; 


N 9 


= E 15 ; 


N w = 


Ei — E 5 ; 


N u 


= E 2 ; 


#12 = 


— E 3 + 2£ 5 


N 13 


= —E 4 ; 


#36 = 


£ 5 . 



(2.19) 

In octet symmetry and at the same order in CHPT, a minimal set of counterterms 
contributing to the AS = 2 component of the K° — K° mixing is: 



(4) 



C Fn G 



AS=2 — ^ A 
+ AA 



27 DxOf 5=2 + D 2 0? s=2 + D 4 Of= 2 + L> 5 0f 5=2 
As ' - 2 + D r O AS = 2 



+ h.c. . 



Notice that this basis is not exactly the one used in 
minimal set. The AS = 2 operators are given by 



(2.20) 
the one in ( |2.20|) is a 



AS- 


=2 


tr (A32X+; 


tr(A 32 x+) ; 




=2 


tr (A32X-; 


tr(A 32 X-) ; 


AS- 


=2 


tr (A 32 u M ) 


tr(A 32 {«",*+}) 


AS- 


=2 


tr (A 32 w M ) 


tr(A 32 K, X _]) ; 




=2 


tr (A32X+: 


tr (A 32 m m u m ) ; 


o AS - 


=2 


tr (A 32 n M ) 


tr(A 32 u^) tr( X+ 



(2.21) 

Since the 27-plet operators with AS = 1 in ( j2.16| ) and the AS" = 2 ones above 
are components of the same irreducible tensor under SU(3)l x SU(3)#, the A 
couplings in both Lagrangians have to be the same. 

The divergences associated with the minimal set of counterterms in ( 2.16| ) 
and (|2.20|) can be extracted from Kambor et al. [11] with the use of the strong 
equation of motion, partial integration and the Cayley-Hamilton theorem. We 
explicitly verified the results in [11] and differ in the 27-plet sector by an overall 
signf||. The subtraction procedure is defined in the usual manner by 



Ei = EZ + 



v 



d-4 



16vr 2 



J^l + \ [lE - 1 - ln(4vr)] 



and 



a = a + 



V 



d-4 



d 



1 1 r 

—4 + 2 [lE 



16tt 2 

In the strong sector we need the counterterms 

z/" 4 f 1 



ln(4vr)] \ 7i . 



16tt 2 



{^4 + 5[*- 1 -M*r)]}r, 



(2.22) 



(2.23) 



(2.24) 



4 EP has independently redone the Generating Functional calculation of the infinities and 
agrees with [pil, |23| and modulo the 27-plet overall sign. 
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with i = 4, 5, 6, 7, 8, and the one associated with the contact term 

with j = 2. The coefficients of the divergent parts are fixed to be 

11 



r - 1 r - 3 r 

1 4 — -, 1 5 — n ) 1 6 



and those of Table |l|. 



144' 



(2.25) 



r 7 = 0, T 8 = — , and 5 2 = — , (2.26) 
1 ' s 48 24 v ; 



Ei 


£i 


e'i 


iV c 


A 


7< 


iV c 


1 


1/4 


5/6 


1 


l 


-1/6 


1 


2 


-13/18 


11/18 


-2/3&L + 0(1/N C ) 


2 





1 


3 








1 


4 


3 


N c 


4 








1 


5 


1 


1 


5 


-5/12 


5/12 


1 


6 


-3/2 


1 


10 


1 


3/4 


N c 


7 


1 


1 


11 


-1/2 





N c 








12 


1/8 





N c 








13 


-7/8 


1/2 


1 








15 


3/4 


-3/4 


1 









Table 1: The divergences and the leading in 1/N C behaviour of the weak 0(p A ) 
octet counterterms Ei and 27-plet counterterms A- Notice that Fq (i.e. an N c 
factor) is factored out. 



2.2 1/N C Counting 

It is also useful to know the 1/N C counting of the different weak couplings in the 
Lagrangians ( 2.16 ) and (|2.20|) . We remind that in this counting F 2 is order N c , 
while i?o is order 1 | 2lfl . The effective operators in (|1.3| ) are four-quark operators 
and the leading contributions to them are of order [0. In particular the 
leading large N c contribution (i.e. in the absence of gluonic corrections) to the 
| AS 1 1 = 1 operator comes from the one W-exchange diagram, while the Box 
diagram leads the leading contribution to the |A5| = 2 operator. In both cases 
only one effective four-quark operator arises in the large N c limit: 



O 



as=i(x) = Q 2 (x) = 4 



{uL7fidl 



(x 



(2.27) 
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1 transitions and 



O 



AS=2 



X 



Q 



AS=2 



X 



4 (s L ^d L ) 



x 



(2.28) 



for \AS\ = 2 transitions. Above we defined qL,(x) = [(1 — 7s)/2] q(x) and sum- 
mation over colour indices is understood inside each bracket. For iV c — > oo both 
currents bosonize independently. So the terms in ( |2.f 6| ) and ( |2.20| ) with two 
flavour traces, each of the current type, may receive contributions of order N 2 . 
However, to assign the correct l/N c counting to single and double flavour traces 
in the weak effective Lagrangian the operators have to be traceless. This can be 
best seen at the quark level. Light four-quark operators are described by SU(3)x 
x SU(3)# tensors 

i^^Tq^T'q^ (2.29) 

with the adequate Dirac structure (pseudo-scalar, scalar, axial-vector and vec- 
tor) 5 . The coupling modulating the traceless part of t, namely t l ^ kl — (2/3) t m ^ ml 
is leading in the 1/N C counting, i.e. of order N 2 . The rest of t is of order N c since 
it has an additional flavour trace which is 1/N C suppressed. Already at order p 2 
the weak octet Lagrangian, if not written in terms of flavour traceless operators, 
does contain double trace terms which are of leading order N 2 . At order p 4 , the 
basis of 27-plet operators in ( |2.18j ), the one in ( [2.21 ) and the one used in Kambor, 



Missimer, and Wyler in JTTJ are written in terms of traceless operators so that 
the 1/N C counting is correct in those cases. 

Neither the octet basis we use in (|2.17|) nor the one in [1]J are written in terms 
of traceless operators so that one has to proceed in two steps to do the correct 
1/N C counting. First, writing the octet Lagrangian in terms of flavour traceless 
operators. The second step eventually needed is the reduction to a minimal basis, 
which involves the use of the eqs. of motion, integration by parts and Cayley- 
Hamilton relations. Again, some of these relations can spoil the correct 1/N C 
counting. In this respect the counting given in [IT] for the octet sector is not 



fully correct, while some of the Cayley-Hamilton relations used in [23, ^4] for the 
reduction to a minimal basis were not appropriate for a correct large N c counting. 

A way of deducing the leading in l/N c contributions to the weak couplings at 
order p 4 is the use of the strong factorization assumption (See [2~3| and references 
therein) where all the weak parameters in ( |2.I6|) are known to all orders in CHPT, 
in terms of the parameters of the strong effective Lagrangian. All the large N c 
contributions to the weak parameters are contained in their factorizable part. 
Up to order p A (and assuming the fudge factor of naive factorization kf — 1) the 
weak couplings receive the following large N c contributions 



G 



27 



5 The tensor i can always be decomposed in a symmetric part ft' — and an antisym- 
metric part ft' kl = — p l - lk . We assume t is symmetric. If t is antisymmetric then t is flavour 
traceless and the couplings modulating it are order N%. 
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4 

Eio = 8L 4 + 2L 5 , En = 8L 4 + 4L 5 , E12 = — 4L 4 — -L 5 , £7 13 = — 4L 4 , 
£> 4 = 4L 5 , D 7 = 8L 4 (2.30) 



and all the others in (|2.16 ) being zero. The strict large N c limit would also 



imply L 4 = 0. The factorization assumption actually corresponds to keeping L 4 
non-zero in (|2.30|) and leaving Gg and G27 as free parameters. 



At next-to- leading order in 1/N C (i.e. in the presence of gluonic corrections) 
all the weak parameters receive corrections both non-factorizable and factorizable 
of short- and long-distance origin. G§ 7^ G27 at this order. Analogously all the 
p 4 couplings in ( |2.30| ) and the rest of the counterterms in ( j2.16|) which are zero 



in the large N c limit do get unknown contributions at next-to-leading order. 
Their estimate is one of the main challenges of low-energy physics. Above we 
have included the factorizable next-to-leading in 1/N C contributions too - the L 4 
parts. As shown in Table [I] the octet counterterms E w , En, E u and the 27-plet 
counterterm D4 get leading order contributions in 1/N C . To arrive at the basis 
we are using in ( |2.16| ) and ( |2.20| ) we made use of the equations of motion, partial 
integrations and Cayley-Hamilton relations. The latter have been used in such a 
way that the 1/N C counting is not broken. The use of the equations of motion 
does not break the counting either. However some integrations by parts do, e.g. 
one can remove the apparent current-current structure of some operators in this 
way. Therefore we have to do the 1/N C counting of the other counterterms in 
Table [l] before the integrations by parts are done. Afterwards, the counting is 
translated to the basis in ( |2.16| )p|. 

We summarize the 1/N C counting of the weak couplings we are using. The 
0{p 2 ) couplings G 8 and G27 are order 1 while G' 8 is order l/N c (notice that Fq 
has been factored out). The 0(p 4 ) couplings E w , En, E12, and _D 4 are order N c 
(Fq has been factored out in (|2.16|) ). The couplings E\, E 3 , E4, E 5 , E13, E± 5 , 
D\, D 2 , D 5 , Dq, and D 7 are order 1. The combination of couplings E 2 + 2Ei/3 
is order 1/N C . 



3 Two-Point Functions 

In this section we give the two-point functions in ( [Lip at order p A for all the 
relevant ij combinations. 



6 In the strict iV c — * 00 limit the singlet 771 degree of freedom needs to be included. This 
one should afterwards be integrated out, leading to counterintuitive 7V c -counting as is the case 
for L7 in the strong sector Since the numerical value of L7 is such that this counting 

seems inappropriate for our real world we neglect this issue. It should however be included in 
modeling approaches since there can be sizeable contributions from it as e.g. seen in Bk @- 
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Figure 1: The diagrams contributing to the AS = two-point functions. A line 
is a pseudoscalar meson propagator, a dot a strong vertex with only meson legs 
and a cross a vertex from the strong Lagrangian with one or more insertions of 
the external pseudoscalar currents. 



3.1 Strangeness Zero 

Here we give the two-point Green's functions in ( |1 . 1|) with % = j to order p 4 . 
Notice that those conserving strangeness with i ^ j vanish since m u = m^. The 
poles of these two-point functions define the masses of the corresponding mesons 
and set the renormalization factors Zj for the pseudoscalar sources P l needed for 
the reduction procedure. 



mt 



(3.1; 



To order p 4 , using m u = = m and neglecting electromagnetic corrections, we 
get from the diagrams in Fig. |l] 



ZjrO = Z n -i 



Zj^o — Zj(+ 



IBlFl 



1 + -2 (2^4 + 774) (4L e 



2/^ - 2/i K - -fM Vs 



2B 2 F 2 



l + -^{2m 2 K + ml) (4L e 



1 

6' 



I? 



2BlFl 



l + -^(2m 2 K + mi) (4L £ 
o 



L 4 ) + ^77C(4L 8 - L 5 ) 

^0 



U) + -^m 2 K (4L 8 - L 5 ) 



L 4 ) + -^ml(4L 8 -L 5 ) 
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2^tt — -/Ak — 2fi va 



(3.2) 



Z 'ir° - Z 'n+ - Z K° - Z K+ - Z ' m ~ 8 ^o( 2 -^8 ~ H 2 



(3.3) 



and 



ml 



m 2 K o 



ni 



</8 



2mBn 



8 



l + ^(2m^ + m^)(2L 6 -L 4 ) 



< + = (m + m.)flb 



l + ^(2m^ + m2)(2L 6 -L 4 ) 



~(m + 2m s )B 



I + (2m 2 K + ml){2L & - U) 



4 



+ ^2 m »?8 ( 2L § - ^5) + - -//„, 



+ 277150 



128 



+ B 2 (m s - ™) (3^7 + L 8 ) . 



(3.4) 



For completeness and later use we also quote the decay constants fx, and 
f m , to the same order: 



fn° — /tt+ 
fx — fx+ 
f 



F 2 



1 + — 2 {2m\ + m 2 )L 4 + —^mlL 5 - 4fi n - 2/i K 
^0 *o 



1 + (2m^ + ml)L A + -^m 2 m L 5 - 6fi K 



(3.5) 



We use the notation 21 



fJLi 



m 2 , / m 2 \ 
In — H . 



32tt 2 F 2 



(3.6) 
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3.2 Strangeness One 



In this section we give the renormalized two-point Green's functions in (|1.1|) up 
to order p 4 . We define them as 



n, 



(3.7) 



Only the analytic contributions from the counterterm Lagrangian Uij(q 2 ) | Count 
are written here, while we give the non-analytic contributions from the one loop 
integration Uij(q 2 ) | Logs in Appendix [A|. The contributing diagrams are shown in 
Fig. |2| With the pseudo-scalar sources P l (x) = P K (x) and Pi(x) = P Vs (x) the 
two-point Green's function is given by 



n^o^ (q 2 



Z K oZ vs 



c 



+ 



{q 2 - m 2 K0 ){q 2 - m 2 m ) y/E 

f 2 4- f 2 
q\G 27 -G 8 + Jvs l lK G' s 



'/s 



F 2 



(2 



j n 



,2 Jrj&_ 

7 o 2 



{2g 4 (G 8 E 3 - G 27 D 2 



, 2 2 

+ q m K 



2E W 
G27 



G 8 (aE x + AE 2 + y E 3 + 8E 4 + 6E 5 
AE U -8E 



12, 



(4D a - 7D 4 -D 6 - 6D 7 ) - 8G' 8 Ul 6 + 7 -L^ 



1 2 2 
+ q 



16 

C; s j 2E 2 - —E 3 - 8E 4 - 6E 5 + 3E U + 8E 12 



G 



27 



(4£>i -D 4 -D 6 + 3D 7 



—m 



K 



G 8 {Ei + £"2 



G 



27 



8Gg (^2L 6 — -Lf 
D x - 8G' (l 6 + h s 



(-Ei — E 2 



G 



27 



3™" 



G 



G 8 E 2 + 2^D l - 8G' [Le-^-L 



D, + 8G' (L 6 — ^-L 



1 



(3.f 



Notice that the renormalized meson masses m n , mx, m m are used everywhere. 
For the case of pseudo-scalar sources P l (x) = P K (x) and P^(x) = P n (x) the 
result is 



U K o v o(q 2 



c 



I Count 



x 



(q 2 - m 2 K „)(q 2 - m 2 ) y/2 
q 2 [ ( 7 



F 4 
2—^ 



r27 



f 2 4- f 2 
U"8 i — — u-« 



rn 2 JjLC' 

-^0 



13 



x ■ X x X 



0_ ^1 



x ■ x x ^—^ X ■ X 



x ■ ^x^ x x 



(*) Q 



¥^ x ^0r x 



Figure 2: The diagrams contributing to the |A,S| = 1 and \AS\ = 2 two-point 
functions up to one loop. In addition to the symbols of Fig. [l|, the full square is 
a weak vertex with only meson legs and the circled cross is a weak vertex with 
one or more insertions of the external pseudoscalar currents. 
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+ - f2 {2q A (G 8 E,-G 27 D 2 ) 

+ q 2 m 2 K [G 8 (4#i + AE 2 - 2E 5 - 2E W ) 

(AD, - 3D 4 -D 6 - 6D 7 ) - 8G' 8 (4L 6 + L 8 ) 



G27 



+ q 2 mi [G s (2E 2 + 2E 5 - E n 
+ 



% (4Di + 3D 4 - D 6 + 3£> 7 ) - (2L 6 + L 8 ) 



2m 2 K ml 



Gs (Ei + £"2) 



—^-Di — 8G' 8 L G 



Finally, for P l (x) = P K+ and Pi(x) = P n ^ we obtain 



G 



G 8 E 2 + 2^fD 1 - 8G' 8 (L 6 + L 8 



(3.9) 



Tl K+w+ (q 2 )\ 



y/Z K +Z^ 



(q 2 - m 2 K+ )(q 2 - m\+ 
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x 2- 



q 2 G 8 



G 



27 



J 77 



F 2 
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+ 
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F 2 



-2q A G 8 E 3 + 2 
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27 



D, 



q 2 m 2 K [G 8 (AE 1 + AE 2 - 2E 5 - 2E 10 ) 

^ (-AD 1 - 2D 4 + D 6 - AD 7 ) - 8G' 8 (4L 6 + L 8 

q 2 ml [G 8 (2E 2 + 2E 5 - E n ) 

% (4£>i - 2D 4 -D 6 - 2D 7 ) - 8G' S (2L 6 + L 8 ) 



+ 2m 2 K ml 



G 8 (Ei + E 2 
G 



G 



27 



D\ — 8G' 8 L 6 



G 8 E 2 + 2^D 1 - 8G' 8 (L e + L 8 )] } 



(3.10) 



3.3 Strangeness Two 

The two-point Green's function for the AS = 2 transition was already calculated 
in |8|]. We include it here for sake of completeness. With the notation used in 
the present work we need the AS* = 2 part of the two-point Green's function in 
(|TID with P l (x) = P K ° and P*{x) = P^°. This gives 

n^-J(g 2 ) (count = — -7—^ — — o x„ Cas=2 G 2 7 4- 



(q 2 -m 2 K0 ) 2 



f 



K° 



q 2 + ^{-M A D 2 



+ 2q 2 (2D 4 m 2 K + D 7 (2m 2 K + m 2 n )) - W 1 (m\- 
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(3.11) 



The non-analytic contributions are in Appendix [A]. There is another contribu- 
tion [AS = ±1) 2 to this two-point function which comes from expanding the 
exponential in ( |1.1| ) up to second order. They are the so-called long-distance 
contributions to K° — K mixing. 



4 K — ► 7T7T Amplitudes 

We have the following decomposition into definite isospin quantum numbers in- 
variant amplitudes [A = —iT], 



A 
A 



K s 
K s - 



71 71 



2 2 

1 A ° + 7S A2 



A 



7T 71 



An 



(4.1) 



Where K s ~ K° + e K° 2 , K° 1{2) = {K° - (+)K°)/ v / 2, and CP K° 1{2) = +(-)#? (2) . 
Since CP violation is small we set e = and therefore QmGg = 0, ^smG 27 = 0, 
and SmGg = 0. We have also included the final state interaction phases into the 
amplitudes A and A%. For the isospin 1/2 amplitude we have 



A = -iao e iS ° 
and for the isospin 3/2 amplitude we have 



A, 



-ML2 e 



(4.2) 
(4.3) 



To order p 2 we get 
a = 



a 2 



^8 + g^27 



V6F {m\ - ml) 
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and 



5 = 5 2 



0. 



(4.4) 



(4.5) 



The order p 4 counterterms contributions to Aq, Af? , and A 2 (see Appendix || for 
the non-analytic contributions) are 
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-cg 8 Vq f ° 4 



fKft 
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77L 
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2 <-AL 4 -L 5 + 8L 6 + 4L 8 ) + 2m 2 Y L 4 
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Fo 2 



(4.6) 



9 iKft 

2ml (-6^1 - 2D 2 + 2D 4 + 6D 6 + D 7 ) 



+ (D 4 - D 5 - 9D 6 + AD 1 



(4.7) 



and 



2 | Count. 



-CG 



27- 



10^3 F 4 



9 f K ft 



1 + 



2m 2 {-2D 2 + 2D 4 + D 7 ) + (D 4 - D 5 + 4L> 7 



(4.t 



The diagrams are depicted in Fig. ^. In addition there are the corrections on the 
external legs and on the internal propagators of the tree level diagrams. 



5 K —> tttt from K —> tt, r]s Amplitudes 

We discuss here the information we can extract from the K — > n,r] 8 two-point 
functions and make some remarks about the parameters needed for K — > tttt 
we cannot obtain. As discussed in Sect. the weak mass term contribution 
to K — ► tttt decays can be absorbed in a redefinition of the other coefficients 
PI , while this is not true in the case of the two-point functions. In typical 



approaches used in lattice QCD or effective models like the one proposed in [ |iq] , 
the weak mass term should be treated as an extra parameter in the determination 
of K — > tttt amplitudes at order p 4 . 

A few other remarks are needed here. The contribution from the 27^ and 
8l cannot be easily disentangled in general, since for m s ^ = m u the two 
components are mixed by higher order effects in m s — m. The A J = 1/2 and 
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-< 









Figure 3: The diagrams contributing to on-shell K — > ttti amplitude up to one 
loop. Symbols as in Figs. [I] and In addition, renormalization of external legs 
and internal propagators of the tree level diagrams have to be added. 



AI = 3/2 contributions are however separate. The two-point function n^o % (g 2 ) 
and the combination (l/v^II/fo^^ 2 ) — n i ^+ 7r +(g 2 ) are pure AI = 1/2, while 
\/2n^o(g 2 ) + Il K + n +(q 2 ) is pure AI = 3/2. 

The K — > 7r and if — > r]g two-point functions defined here are not measurable 
in experiments, so obtaining the counterterms from them implies that one has 
to calculate the relevant two-point functions either using lattice QCD or 
using other hadronic approaches || || . 

Recent work on AI = 3/2 transitions in quenched CHPT \27\ gives numeri- 
cally consistent results with quenched results on Bk- The main uncertainty are 
the unknown quenched counterterms: G27 and D^s. The value of several of them 
can be similarly extracted from two-point functions but this has not been done 
so far f28|. In the near future it might be possible to calculate the two-point 



functions we propose in ( |1.1| ) both in the quenched and unquenched case and the 
K — > Ten amplitude only quenched. The comparison of all the constants that can 
be calculated using the simpler two-point correlator would then be a useful tool 
to estimate quenching errors on the remainder. 

We now discuss the expressions of Sect. ^| to check which constants are ob- 
tainable. The order p 4 counterterms Ei 3 , E\§, and D 5 cannot be obtained from 
the Uij(q 2 ) since they do not contribute to them. They do however contribute to 
K —>■ 7T7T. Their value in the large N c limit and with the factorization assumption 
is known, see Sect. |2~2 . 

In the chiral limit we can get G§, G27, E 3 , and D 2 . Away from the chiral limit, 
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we can get G' 8 from the terms quadratic in the meson masses. From the terms 
quartic in the meson masses we can get E\, E 2 , and D\. The terms proportional 
to q 2 m 2 T j < allow us to obtain the combination 2E W + En, _D 4 , _D 6 , and D 7 . The 
latter also determine two more combinations of couplings, though not needed for 
K — > 7T7r: £4 + E12 and E 5 + E w . E w is needed for K — ► 7T7t but cannot be 
separately disentangled from the two-point functions Uij(q 2 ). 

We have determined three order p 2 couplings and eleven order p A ones. We 
have in addition four relations which are independent of the value of the couplings 
to test chiral symmetry at order p 4 . This is the main result of this manuscript. 
We are left with four unknowns E w , Ei 3 , Ei 5 , and D 5 in order to extrapolate to 

K — > 7T7T. 

What can be said about the missing coefficients. Three kind of arguments 
can be used: 

1) Order of magnitude: We know the leading in 1/N C contributions to all of 
them, Eq. ( |2.30| ). In particular, Ei 3 , Ei 5 , and D 5 only receive non-factorizable 
contributions. From the discussion above, we have seen that we can determine 
eleven next-to-leading in 1/N C contributions to the couplings, this should give us 
some information on the ones we cannot get. 

For instance, assuming that all the 1/N C contributions are of the same order, 
since 2E W + En is obtained from terms of the type q 2 m 2 K and we know its 
factorizable contribution 

{2E\q + En) 

Factorizable 8L 5 + 24L 4 , (5.1) 

we can use (conservative choice) the next-to-leading result we have for 2E W + En 
as a good estimate for Ei 3 , J5 15 and the non-factorizable part of E w . Analogously 
for D 5 we can use the result 

(AO Factorizable Factorizable 0J-/4. , (5.2) 

for predicting the non-factorizable part of D§. 

2) Resonance Saturation: Another possibility is using estimates of higher or- 



der parameters coming from resonance exchange saturation, as done in [p4[ . It 
is well known and experimentally proven that vector and/or axial- vector domi- 
nance is not at work in the weak sector (which is instead the case in the strong 
one). Within the tensor formulation of vector (axial- vector) resonances used in 
[ p4}1 all the counterterms in ( |2.16| ) only receive contributions from scalar and/or 
pseudo-scalar resonances. In Table § we summarize the resonance contributions 
to the counterterms contained in ( |2.16| ). The notation is the one used in 



For the 27-plet the derivation is done in Appendix 0. All of the unknown terms 
Ei , Ei 3 , Ei 5 and D 5 only receive scalar and/or pseudo-scalar contributions. This 
is consistent with the consequences of the factorization assumption as shown in 
(2.30), where the weak counterterms are all expressed in terms of the strong 
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c m9s 
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Table 2: The contributions to the octet Ei and 27-plet Di counterterms in ( |2.16| ) 
from scalar octet (S), scalar singlet (Si) and pseudo-scalar octet (P) resonance. 
The pseudo-scalar singlet (Pi) only contributes to E4 with a term d m gp. A 
factor 1/M^ is pulled out. Vector and axial- vector resonances exchange do not 
contribute p4|. 



counterterms L4 and L5, which are in turn saturated by scalar exchange |24|, 
A simplified resonance model is e.g. the one where pseudo-scalar resonances ex- 
change is neglected []. Their contribution is small in the strong sector [29 . With 
this reduction several relations are valid. E s = 0, E 4 = 0, D 2 = and D 5 = are 
a test of scalar dominance, while 



2E[ + E n 



— (2E[ 

Cm 



1E\ + 3Ey 







(5.3) 



7 One can, of course, use other resonance models to make the analysis. We mention that 
in the vector formulation with vector fields used in pc[ , not antisymmetric tensor fields as in 
]24| , i?i5 also receives contribution from vector resonance exchange. We do not address here 
the question of the equivalence of different resonance models in the weak sector. 
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tests the absence of singlet scalar resonance contributions. The couplings Cd and 
c m are scalar couplings from the strong sector [24, |29j. Additional relations are 
(c m /cd) E[ 5 + El — El = 0, valid in the absence of pseudo-scalar resonance 
contributions, and 



E 



1.3 



-(2E r w + E r xl ) 



c 



(E{ - El) + — (3JVJ + 3iV 2 r + 6jV 3 r ) 



(5.4) 



valid in the absence of the singlet scalar resonance. The combination of TVj 
corresponds to K% which receives no vector /axial- vector contributions and is a 
combination of the p 4 constants that can be determined from K — > 3ir decays 
p4| , Hence, introducing additional information about K — > 3ir decays we 
can use resonance arguments to determine Ei 3 . One main observation is that the 
direct determination of most of the couplings from the analysis of the two-point 
functions offers already a powerful test of the validity of different resonance sat- 
uration assumptions. 

3) Factorization: We can of course also adopt the strong factorization assump- 
tion as often used to get at the undetermined parameters. Again this procedure 
can be well tested by the fourteen parameters we can actually determine and 
comparing them with the predictions of Eq. ( |2.30| ) with G 8 and G27 as free 
parameters. 



6 The Weak Mass Term Contributions 

In the literature there are conflicting opinions about whether the weak mass term 
contributes to K — > 7r and K — > vacuum matrix elements. In jl5] the claim is 
that they do not and in @, |7| they do. The underlying reason for the difference 
is that the meson fields used in those two references differ by a field redefinition. 
For on-shell matrix elements this makes no difference and as a consequence both 
analyses agree for the K — > tttt amplitudes. Since neither the K — ► 7r transition 
nor the K —>■ vacuum transition can be allowed on-shell, if the masses are such 
that K — > tut is possible on-shell, we first have to correctly define what we mean 
by off-shell matrix elements. 

Green's functions defined by quark currents, as introduced in CHPT by Gasser 



and Leutwyler [21], are well defined for all values of momenta and thus provide a 
proper definition of off-shell quantities. We have shown here that for the pseudo- 
scalar currents as sources the two-point function that defines properly an off-shell 
K — > 7r transition does depend on the weak mass term, or the coefficient G' 8 . 

Of course, we also find that for the on-shell transition K — > tttt, the weak 
mass term does not depend on G' 8 to order p 2 as shown before in 0, 14, 15]; but 
it does have contributions at higher order, see the discussion below. 



The discussion of |11| and 12(1 indicates at which level the weak mass term 



can contribute to K — > txtx amplitudes. The basic argument can be phrased in 
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terms of the strong equations of motion for the u field in ( |2.2| ). Terms that can be 
removed using the equations of motion can be described by changes in the other 
coefficients of the Lagrangian. The argument of Sonoda and Georgi is that the 
weak mass term can be written as a total derivative assuming m s 7^ rrid (see also 
[pH[] ) and as such does not contribute to physical amplitudes. This argument fails 
in the presence of external fields. It crucially requires s + ip = A4, otherwise other 
non- derivative terms involving s + ip — M. remain after the use of the equations 
of motion. That is why the argument fails for the non-tadpole diagrams like the 
two-point functions considered here. These extra non-derivative terms give the 
G' 8 contributions to our ILj(g 2 ) two-point functions. 

To see what happens in the case we add external vector and axial vector 
fields, let us show the argument of Sonoda and Georgi in more detail extended to 
include vector and axial-vector external fields. The basic underlying argument is 
that the effects of terms that vanish using the lowest order equation of motion 
can be described by changes in the other parameters of the effective Lagrangian. 
As we will show below this implies that the effects of the weak mass term can 
be absorbed in shifts of the other parameters for all processes involving on-shell 
pseudo-scalars and photons. As a consequence the contributions from the G' 8 
term vanish at order p 2 for these type of diagrams and can be absorbed in shifts 



of the other parameters at higher order, i.e. the Ei and Di of |TT| at order p 4 . As 
stressed earlier, this does not mean that the contributions from the weak mass 
term are zero in this case, only that they can be described by shifts of the other 
parameters. 



The equation of motion from the lowest order Lagrangian in (|2.2|) is 



2D^D»U) - tf X + X ]U ~\ ^(-U j X + X ] U) = . (6.1) 
When the external scalar and pseudoscalar fields are zero this becomes for i ^ j 

2D^D il U) ij -2B^j\ 3 m 3 ^2B Q U ij m i = 
-2D ii (JJD^) i ,-2B^j\ 3 m l ^2B {) U l jm ] = 0. (6.2) 

The second equation can be derived by first multiplying Eq. (|6.1|) on the left by 
U and on the right by W and using the unitarity of U extensively. The weak 
mass term for AS = 1 transitions is proportional to 



( X f U + C/t x ) 23 = 2B {m d U 23 + m s Ul 
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= 2^±^(D^D^U)) 23 + 4J%^iDp(UD*Ul)) m . 
nig — m d trig — m d 

(6.3) 

For external vector and axial-vector fields zero or equal to the photon field the 
last line is a total derivative and thus does not contribute to the action. This 
proves the comments made above. 
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Figure 4: The Feynman diagrams for K$ — > 7r°77 at order p 4 . Notation as in 
previous figures. The strong vertices can now also be the Wess-Zumino term. 

Processes with on-shell pions and kaons and photons could in principle depend 
on Gg already at their lowest non-zero order. We have verified that G' 8 does 
contribute to K L — > 77, and hence to K L — > fi + fi~, at lowest non-zero order (i.e. 
at order p 6 ) similarly to the part from GgQ. For Kg — > 7r°77 there is already 
a contribution at order p 4 . The relevant diagrams are depicted in Fig. |j. The 
contributions from the weak mass term are non-zero but local, the non-local 
parts containing — fn 2 M ) cancel, and thus as proved in general above can 

be absorbed in shifts of the Ei. 

These results are important since this contribution has never been included 
in the long distance estimates of the processes above. Notice also that the value 
of G' 8 is unknown. 

In addition, physical amplitudes like K — > tctt that do not depend on G' s at 
lowest order can depend on it at next-to-leading order in CHPT. In this case 
the G' 8 contribution can be reabsorbed in shifts of order p 4 and higher couplings. 
Nevertheless notice that in analytical predictions, the value of G' s should be known 
to do this shift. 

7 Conclusions 

We have extended the p 2 analysis of [|J with the non-analytic contributions of 
H in two ways: first, we have included all the relevant p 4 couplings and the full 
one loop contributions to the two-point functions K —> n,r] 8 and the K —* nn 
amplitude and second, we have changed from the vague notion of an off-shell 
meson field to the well-defined notion of Green's functions of external fields in 
the presence of the weak non-leptonic interaction. 

We have confirmed the results of [Q that the three relevant couplings at order 
p 2 can be fully determined and discussed the necessity of including the weak mass 
term in this analysis. 

We concluded that the weak mass term can contribute in some physical am- 
plitudes even at their lowest non-zero order and have given several examples 

8 See |}2| for the standard discussion and the warning about SU(3) breaking parameters 
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where this happens. In general this can happen in the case of on-shell Green's 
functions which receive contributions from off-shell flavour changing two-point 
functions. As an example, the weak mass term gives lowest non-zero order un- 
known long-distance contributions to processes like Kl —> 77, K$ — > 7r°77 and 

To order p A in the chiral expansion we find six more parameters in the 27-plet 
sector that could in principle contribute to K — > tt,i] s ,tttt. Of these, we can 
directly determine five from the two-point functions, leaving one, D 7 as a free 
parameter. This parameter vanishes in the large N c limit and is proportional to 
L4 in the factorization model. If the predictions of this model turn out to be 
satisfied by the other five parameters we can take the factorization prediction for 
D 7 and obtain a value for K — > tttt. 

In the octet sector, there are ten more operators of which we can also directly 
determine six combinations. We can use them to test the predictions of various 
models like factorization, the weak deformation model |24j], resonance models, etc. 



The weaker assumptions of resonance saturation by vector, axial- vector and scalar 
resonance exchange allows to determine one more combination of counterterms 
from the two-point functions and one more parameter (-E13) can be fixed if one 
of the slope parameters of K — ■> 3tt, namely K2 of Ref. fl3|| , is known. To obtain 
the full set of counterterms at order p 4 we need to use factorization or another 
more restrictive model. Factorization and alternative models can be strongly 
constrained by the value of the parameter combinations that can be directly 
extracted from K — > 7r,i] 8 two-point functions. 

As a calculational tool, we have provided the complete one-loop formulas for 
the two-point functions of the octet symmetry case and for the on-shell K — > tttt 
amplitudes with quark masses m s 7^ = m u all different from zero. In |13[] and 



[[24[] the combinations suppressed by m^/m 2 K were neglected. This might be a 
good approximation for the real quark mass values [at the level of 10 % though], 
but will not necessarily be true on the lattice. 

Since our two-point functions are much easier to determine on the lattice they 
also provide a better laboratory to study unquenching effects in the non-leptonic 
weak sector than the full K — > tttt amplitudes. 
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B Order p A Loop Contributions to K —> tttt Am- 
plitudes 

In addition to the definitions used above we need now, 
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C The 27-plet Weak Lagrangian from Resonance 
Exchange Saturation 

We derive the weak effective Lagrangian at order p 4 in the 27-plet sector for 
AS = ±1 transitions by assuming resonance exchange saturation of the couplings. 
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We restrict the derivation to those terms listed in ( |2.18| ). The derivation of the 
octet sector can be found in [24}|. We refer the reader to Refs. |29| for 



details on the method. The weak 27-plet Lagrangian in ( 2.16Q can only receive 



contributions from scalar (octet 5* and singlet Si) and pseudo-scalar octet P 
resonances. The relevant weak couplings of the light meson fields to resonances 
can be written as follows: 

^ 7 = E^f +Y.?pK( + f s K^ +h.c, (C.l) 



i=i 



where 

Kf = 4 tr(A^) tr(A w x+) , = t% tr(A i3 -{S, tr(A fc ^) , 

K$ = t% tr(A^) tr(A«uX) > 

K[ = it% tr(AyP) tr(A fc/X -) , = it% tr(A^, P}) tr(A w u") , 

K Sl = Si tr(A^) tr(A fe ^) . (C.2) 

Inserting the lowest order solution of the equations of motion for the resonance 
fields in ( |C2| ), the 27-plet weak effective Lagrangian at order p 4 and order Gf is 
given by 

4?= E (C.3) 

R=S,P,S 1 1V1 R 

where J,f , j£ are the strong and weak currents respectively, coupled to the res- 
onance R at lowest chiral order p 2 . The weak currents are defined from ( |C.1|) as 
follows 

C% = tt(SJi) + tr(Pj£) + Sis** , (C.4) 
while the strong currents are given by [^9| 

J* = °d + c m X+ Js 1 = Cd tr(u^) + c m tr(x+) 

Js = idm X- ■ (C.5) 

The contributions to the low energy 27-plet weak effective Lagrangian (|C.3|) are 
summarized in Table 0. 
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